Nonmonotonic inelastic tunneling spectra due to emission of surface magnons in 

itinerant electron ferromagnetic junctions 
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We propose an itinerant electron theory of nonmonotonic inelastic tunneling spectra of ferromag- 
netic junctions observed in several recent experiments. Assuming that the Coulomb interaction 
between electrons across the barrier is poorely screened due to the exponential suppression of the 
band electron density inside the barrier, we consider a process where a hot electron from one ferro- 
magnet during the tunneling excites a magnon at the surface of the other ferromagnet in the course 
of the exchange scattering with a Fermi sea electron from the opposite side. For this process one 
can derive an inelastic transfer Hamiltonian in which the matrix elements of the electron-magnon 
coupling depend on the magnon wave-vector (and hence on its energy) via the long-range Coulomb 
potential. This dependence can account for the sharp peaks observed in the low-temperature tun- 
neling spectra of a number of ferromagnetic junctions at the excitation energies much smaller than 
the Curie temperature of the ferromagnets. 

PACS numbers: 72.25.Ba, 72.25.Rb, 73.40.Gk, 75.30.Ds, 75.30.Et 



I. INTRODUCTION 

Spin polarized transport in tunnel ferromagnetic junc- 
tions has been a subject of intense research"'^ motivated by 
the desire to develop a form of electronics which utilizes 
the dependence of the junction resistance on the spin- 
polarization of carriers in the ferromagnets. Since fer- 
romagnetic metals have more band electrons of one spin 
polarization (known as majority carriers) present at the 
Fermi energy Ep than of the inverse polarization (mi- 
nority carriers), the resistance depends on the relative 
orientation of the magnetic moments in the ferromag- 
nets which is controlled by an external weak magnetic 
field. With parallel magnetizations, the tunneling oc- 
curs between majority (and minority) bands whereas in a 
junction with antiparallel magnetizations carriers tunnel 
from majority to minority bands (and vice versa). The 
resulting spin current mismatch produces a larger con- 
tact resistance in the antiparallel case, an effect known 
as junction magnetoresistance (JMR)2i^. 

This simple picture was successful in interpreting 
current-voltage (I-V) characteristics at low bias voltage 
of the order of a few millivolts. For V about a hun- 
dred millivolts the JMR was found to be strongly re- 
duced^i^i^i^. Since this energy scale is of the order of 
that of spin excitations (about the Curie temperature Tc 
of the ferromagnets), it was suggested that electrons tun- 
neling from one ferromagnet with energies ~ \sV\ above 
the Fermi energy in the other ferromagnet (hot electrons) 
excite both surface^, and bulk?- magnons, so that the re- 
duction in the resistance for the antiparallel alignment 
turns out to be bigger than for the parallel case lead- 
ing to a decrease in the JMR. Within the framework of 
the transfer Hamiltonian model"'-*''^'': and neglecting the 
energy dependence of the electron-magnon coupling, one 
can identify the nonlinear magnon-assisted contribution 
to the tunneling current by taking the second derivative 
d^I/dV^ which can be related at zero temperature to 



the magnon density of states O at the excitation energy 
\eV^: 



d^I/dV^ oc sign(\/)r2(|ey|). 



(1) 



where e is the electron charge. 

In the case of the bulk magnons, Q cx \eV\^^^, the 
inelastic contribution vanishes for V ^ whereas 
for the surface ones the density of states is energy- 
independent, fl = const, and the second derivative 
is discontinuous at = 0. This corresponds to a cusp- 
like zero-bias anomaly in the conductance observed at low 
temperatures in a number of experiments^i^iSii^ pointing 
to the dominant role of the surface magnons in the inelas- 
tic tunneling spectroscopy at excitation energies smaller 
than the Curie temperature. 

In contrast with simple formula Q real inelastic tun- 
neling spectra are strongly nonmonotonic*. According to 
the data obtained in Ref. 8 for Co/Al203/Ni8oFe2o junc- 
tions, d^I{V)/dV^ has two antisymmetric sharp peaks 
at excitation energies much smaller than the Curie tem- 
perature of the ferromagnets. Similar results have been 
recently obtained in Refs. 12. 1,1^1^ for several types of 
junctions, showing that the nonmonotonic behaviour of 
the tunneling spectra is a generic feature of inelastic 
tunneling between ferromagnets. In this paper we pro- 
pose a microscopic theory of magnon-assisted tunneling 
which links the nonmonotonic tunneling spectra to the 
energy dependence of the electron-magnon interaction. 
The magnon emission (absorption) processes are treated 
in th e rnode l of band electrons in the similar spirit as in 
Refs. ll5llfil 

To anticipate, as the excitation energies of interest are 
quite small, we focus on the exchange processes accom- 
panied by emission of surface magnons. They are snown 
schematically in Fig. 1 for the simplest case of fully- 
polarized ferromagnets with antiparallel magnetizations 
where elastic tunneling is forbidden. We assume that, 
due to the exponential suppression of the electron den- 
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FIG. 1: Exchange-induced spin-flip tunneling between half- 
metallic ferromagnets with antiparallel magnetizations: (a) 
- while tunneling, a hot spin-up electron from the left-hand 
ferromagnet excites a Fermi sea spin-down electron from the 
opposite side into a state above the Fermi level. Although 
the tunneling electron cannot occupy any one-particle state 
in the right-hand ferromagnet, its presence on the right-hand 
side can be treated as an excitation of the collective spin 
of the Fermi sea (spin-down) electrons, (b) - temperature- 
stimulated counter exchange process due to the finite occupa- 
tion number of the initial state with energy e + \e\V + uiq^^ 
above the Fermi level in the right-hand system, qy and 
Wqy < \e\V are the wave-vector and energy of a surface spin- 
wave; a and U are the barrier thickness and height, respec- 
tively. 



sity inside the barrier, band electrons interact directly 
across the barrier via a long-range Coulomb potential 
which depends only on the distance between electrons in 
the interfacial plane with the screening radius rs greater 
than the barrier thickness a. In the process shown in 
Fig. la the magnon emission occurs in the course of 
the exchange scattering between a hot spin-up electron 
from the left-hand ferromagnet and a Fermi sea spin- 
down electron from the right-hand system. The latter 
is excited into a state above the Fermi energy, leaving 
a hole inside the Fermi sea, whereas the former cannot 
occupy any one-particle state in the right-hand system 
because they are all spin-down ones. Collectively, how- 
ever, the state of the "spin-down" Fermi sea with an 
electron having the wrong spin can be described as a 
spin-wave excitation^- with wave-vector qy parallel to 
the boundary and energy < \e\V . The contribution 
of such processes into the inelastic tunneling spectrum 



FIG. 2: Inelastic electron tunneling spectrum at different tem- 



peratures: (A) ksT/Dry = 0, (B) kBT/Drl 



0.1, (C) 



ksT/Dr-^^ = 0.3, (D) ksT/Drg^ = 1, (E) ksT/Drg' 
where Dr^^ is the energy of a surface magnon with the wave- 
length equal to the screening radius. 



is characterized by the matrix elements of the electron- 
magnon coupling which depends on the magnon wave- 
vector via the Fourier transform of the Coulomb poten- 
tial V(q||) = 27reV(r52 + 52)i/2^ 



T.V^O 



v^(qii) 

^ V2(0) 



-1 \ " 



de X 



(2) 



X [n'ie)n'ie + \e\V ~ u;^^^) - n'{e + \e\V + c^q„ )n'{e)] . 

In the first term the Fermi occupation numbers n(e) and 
n(e + \e\V — ij^qn) correspond to the initial (i) and fi- 
nal (f) electron states in the process shown in Fig. la. 
The second term comes from a similar exchange process 
(Fig. lb) where in the initial state a spin-down electron 
has energy e + |e|V^ -I- cuq^^ above the Fermi level in the 
right-hand system whereas in the final state a spin-up 
electron has energy e with respect to the Fermi level in 
the left-hand system. For positive bias voltage V > this 
process is only possible at finite temperature T when the 
occupation number of the initial state n{e+ |e|y-|-a;q|| ) is 
not zero. In Eq. Q = A/AttD is the density of states 
of the surface magnons (proportional to the contact area 
A) and ujq^^ — Z3g| is their dispersion relation with D 
being the spin stiffness. With magnon wave- vector de- 
pendence of the electron-magnon interaction, the second 
derivative of the tunneling current Q has two antisym- 
metric peaks as a function of bias voltage (Fig. 2). As the 
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temperature increases, they broaden and shift to higher 
energies in a similar way as in the experimenta^*Sii2ii^ill. 
The peak intensity in the hmit of zero bias and temper- 
ature, \(Pl/dV'^\T,v->o can be expressed in terms of the 
barrier transparency and electron band parameters of the 
ferromagnets. 

The paper is organized as follows. In sectionslTlland llTTI 
we propose a new simple method of deriving both tunnel- 
ing current operator and transfer Hamiltonian for inter- 
acting electrons. Unlike the previous approaches to this 
problemi2ii2iS£, we derive effective microscopic boundary 
conditions for the "right" and "left" electrons at the bar- 
rier walls and show that they can be rewritten in the form 
of the transfer Hamiltonian in the coordinate representa- 
tion. In section IIVI the results for the tunneling current 
between ferromagnets of arbitrary spin-polarization are 
discussed for both antiparallel and parallel alignment of 
their magnetic moments. 



II. EFFECTIVE BOUNDARY CONDITIONS 
FOR THE "LEFT" AND "RIGHT" ELECTRONS 

In this section we find an approximate solution of the 
equation of motion for the field operators of the interact- 
ing electrons inside the barrier and eliminate this region 
by deriving effective boundary conditions for the left and 
right operators at the barrier boundaries. 

Consider a contact of a large area A between two fer- 
romagnetic metals separated by a tunnel barrier charac- 
terized by the thickness a and the length of the electron 
penetration A — h/{2m.Uy/^ which depends on the bar- 
rier height U measured with respect to the chemical po- 
tential and the electron effective mass m (Fig. 1). For U 
of the order of the Fermi energy Ep in the metalsi^i^, one 
can neglect the energy and momentum dependence of the 
electron penetration length. The thickness of the barrier 
a is normally much greater than A and the density of 
band electrons inside the barrier, c^'^^^hq is exponen- 
tially reduced compared to that in the leads no, leading 
to direct electron-electron interactions across the barrier 
via a long-range Coulomb potential. Unlike the case of 
Refs.TE'lff, here we deal with tunneling accompanied by 
the interaction inside the barrier. 

In order to treat such inelastic processes, we will use 
the Coulomb potential V(R) = e^Q-^/^'^/R with the 



Hamiltonian 



screening radius i 
that in the metals 



l/(e "/^rto)^^^ much greater than 
), where R = {x; r} with 
the X-axis perpendicular to the interface and the posi- 
tion vector in the interfacial plane r. In the case of the 
experiments described in Refs. 00, rg can be estimated 
as being greater than the thickness of the barrier a, which 
allows one to neglect the x-dependence of the Coulomb 
potential: V(r) = e^e^''/'"=/r. Out all of the inelastic 
processes we will only take into account those due to the 
exchange scattering of two electrons in the states with op- 
posite spins a and —a since they are known to result in 
the exchange-induced spin excitations"'^-^. The exchange 



^^^^^^ I dx,dx2 X (3) 



Hex = — 

a,k27^ki+q 

2 

Xxlk2-q(^l)X-«k2(a;2)xLaki+q(2;2)Xaki(a;i), 

is written in terms of the operators Xak.{x) and x^i(.(a;) 
annihilating and creating, respectively, an electron with 
spin a and wave- vector k parallel to the interface at point 
X inside the barrier (the index || is dropped from now on). 
The Fourier transform of the Coulomb potential is given 
by 



Vq = 2Tie^ I rdr Joirq) 



27re^ 



irf + q^y/^' 



(4) 



where Jo{x) is the Bessel function. The equation of mo- 
tion for Xak{x) can be written as 



[9^-A-2]x.k(x)=-|^5:vk-k. 



(5) 



kiq 



X / dxiX-ak+q(a;i)xL„ki+q(2;i)Xaki(a;), 



with usual boundary conditions imposed by the continu- 
ity of the particle current at the barrier walls: 

Xak(±a/2) - 9,Xak(±a/2) = d^^^^^l, (6) 

where ^I*^'^ = 4'Q.k(±a/2) are the operators of the right 
and left systems acting at the barrier boundaries. The 
eigenstates in the left and right systems can at this stage 
be arbitrary. 

To first order in the interaction the solution of Eq. 

is 



X«k(x)=xi"^(x)-E^ / 

kiq a 



dxidx2 X (7) 



Xt;(°)(^,^l)xLVq(^2)xTk,+q(^2)xL"^(^l), 

where the operators of the non-interactive system 



xL°,^(2:) = 0'^(^)*^.k + '/''(^)* 



iki 



sinh(a/2±a;)/A 
sinh(a/A) 



(8) 
(9) 



are linear combinations of the two fundamental solu- 
tions of unperturbed equation jSJ. Q^^^ is the Green 
function of the unperturbed equation which can be con- 
structed using (x) as follows: 

gio) 

{xi,X2) ~ —2rnh ^Asinh(a/A)x (10) 

(/)'(a;i)(^''(x2) xi > X2, 
0'(a;2)0''(xi) X2 > xi. 
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Note that our choice of the integration constants ensures 
that solution ^ matches the operators of the right 5"^^ 
and left VP^j^. systems at the boundaries x — a/2 and 
X — —a/2 (first boundary condition in Eq. ©)• 

Solution lO is expressed in terms of 'f j^i^ and ^'Jjj. and 
contains no more constants to be determined. Inserting it 
into the boundary conditions for the derivatives and 
evaluating the integrals over the coordinates, we find 



a iTr'" _ ok ^ ak 



kiq 



TVk-ki 

Ah 



*-ak+q*-aki+q ' * -ak+q * -aki +q J ^aki'K^^J 



rt 



k I „-a/A TVk-ki 

Ah 



kiq 



f*' , , ^'t +^''~ , , W"", (12) 

\^^-Qk+q^-aki+q ' ^-ak+q^-aki+q^ ^akii 



T = maX/h = a{m/2Uy^^. 



(13) 



These equations now serve as effective microscopic 
boundary conditions for the right and left operators. Be- 
cause of the tunneling, the right and left operators are 
mixed in Eqs. (dJ and In the limit — > 0, the 
coupling vanishes: 



As in the previous section, we will not use any particular 
eigenstates in the left and right systems. 

As the boundary conditions (|ll|l and \i2\ conserve 
the current density, the total current operator / can 
be related to the current density operator at any of 
the boundaries, for example, at the left one: / = 
i|Eak(5-*ak*Lk - h.c.) with the derivative d^^'l^ 
given by boundary condition p2|) . It can be written 
as the sum of an elastic and an inelastic contributions: 
I ^ lei + I in , where 



f leT s-^,^r\ 1 , N ^ h^e~°-/^ , , 

= -rE(^I*^k-/^-c.), r=— (15) 



~ _ ieT y-^ TVk2-ki-q 

- h ^ Ah 

Q!,k2#ki+q 



(16) 



2 \^ Qk2-q ctk2 — ctki+q ctki ' 
^Qki-q^-Qki ^-ak2+q^ak2 



The product ^^^k^ak ™ ^'i- represents the energy 
operator associated with a one-particle elastic "hopping" 
between the left and the right systems. 

Note that for electrons with close enough wave- vectors 



ok 



A 



(5 * , = - 



ak 



A 



(14) 



and the boundary conditions describe isolated right 
and left systems with the particle current vanishing at 
both X — a/2 and x — —a/2. Boundary conditions 
(|14|l will be used later to introduce the eigenstates in 
the isolated right and left systems. When evaluating the 
integrals over the coordinates we have only taken into 
account terms linear in e~°'^^ ^ 1 in Eqs. (|ll(l and (|12(l . 

Note that the exchange interaction results in the mix- 
ing of the left and right operators with opposite spins in 
the boundary conditions 111|) and (|12|l . which takes into 
account inelastic spin-flip processes during the tunnel- 
ing. These terms are proportional to the tunneling time 
T assumed to be sufficiently short to justify the use 
of perturbation theory. 



III. TUNNELING CURRENT OPERATOR AND 
TRANSFER HAMILTONIAN FOR 
INTERACTING ELECTRONS 



|k2-ki| «(r^^ + g^) 



2N1/2 



(17) 



the matrix elements of the interaction in Eq. are 
independent of both k2 and ki: Vka-ki-q ~ Vq. There- 
fore, the sum over k2 picks up the products of the oper- 
ators describing simultaneous creation of a hole and an 
electron with opposite spins. The superpositions of such 
electron-hole pair operators can be related to collective 
spin operators of itinerant electrons penetrating into the 
barrier. Let us introduce first the operator of the to- 
tal spin of the electrons penetrating into a thick barrier 
(a ^ A), say, from the right ferromagnet: 



a/2 



2^2 

k 



(18) 



Here we use the effective boundary conditions 1)11(1 and 
(|12|l to derive the microscopic tunneling current operator 
and transfer Hamiltonian for interacting electrons. Com- 
pared to the method pioneered by Prange^^, in our ap- 
proach both current operator and transfer Hamiltonian 
will be expressed in terms of the field operators taken at 
the left, 4"^^ and right, ^'J^^^. boundaries of the barrier. 



Here X"^^{x) JSJl exponentially decays over the distances 
of the order of A from the right boundary. At the end of 
this paper (see, Eq. I0lj) we calculate SI and show that 
it is a macroscopic quantity (much greater than unity), 
so that it can be treated as a classical spin. Then the 
operators S^~^ and , raising and lowering the total 
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spin S^, can be introduced as 

a/2 



SI 



^ 

a/2 

(^:^E^<-A- (20) 



They are normalized in the usual way to satisfy the boson 
commutation relation: 5'q+5!l~ — 5!lq5'q+ — sign{Sl). 
Thus, and S^S^ are the magnon annihilation and 
creation operators, respectively, for positive and vice 
versa for negative S^. 

Keeping in equation only the coherent terms with 
wave- vectors satisfying inequality (|17|l . one can express 
the inelastic current in terms of the magnon operators 
(PJ and HOI) as 



h ^ Ah 

kq 



h.c. 



ik+q*Tk^ 

(21) 



where for simplicity \Sl\ = \S[} = \Sz\- Since the op- 
erators S'q'"'' and 5'q''~ change the electron spin in the 
layers of the thickness A near the barrier boundaries, the 
magnons in question are surface ones with wave-vectors 
parallel to the interface. In what follows we adopt a har- 
monic time evolution of the magnon operators e^^'^'i 
with the parabolic dispersion ojq' = D^'^'q^, where is 
the spin stiffness. 

In principle, one can solve the equations of motion in 
the electrodes with effective boundary conditions Hll() 
and p2(l treating the tunneling coupling in Eqs. 1111) 
and p2|l as a perturbation. Using then these solutions, 
one can average the products of the left and right op- 
erators in the equations for the tunneling current H15() 
and (PT|l . However, instead of this tedious procedure, it 
is convenient to derive an effective transfer Hamiltonian 
and use it for the averaging. 

As the tunneling current operator / is proportional to 
the rate of change of the particle number operator of one 
of the systems (e.g the left one), the transfer Hamiltonian 
Ht can be formally associated with the operator / in the 
standard way: 



i=-eNL = 'UNl^Ht], 
Nl = Eak / <k(2^)*c.k(a;)dx. 

2:<-a/2 



(22) 



One can easily verify that if the transfer Hamiltonian is 
taken in the form 



(23) 



ak 



kq 

J^-^lT,■''t 



the commutator [N^ , Ht] results in the current operator 
/ — lei +Iin given by Eqs. ((TB|l and l(5lj. A more tedious 
calculation allows one to derive Ht directly from the ef- 
fective boundary conditions l(TT|l and l(T^ independently 
of the current operator / and make sure that Ht and / 
are automatically related by Eq. (|22() . which is simply 
the consequence of the current conservation expressed by 
the boundary conditions l|ll|l and l|12|l . 

The first term in Eq. (|23|l is the coordinate version of 
the well-known elastic tunneling Hamiltonian of Ref. 
Expanding the operators 5*^'^ = 5'Qk(±a/2) in the eigen- 
states of the isolated right and left systems, one can go 
over to the momentum representation used in Ref. Hol 
The advantage of the coordinate representation is 
that the tunneling matrix element T pSIl is a constant, 
which makes perturbation theory in the coordinate rep- 
resentation simpler. 

The second term in Eq. H23|) describes inelastic elec- 
tron tunneling accompanied by the magnon emission (ab- 
sorption). It is physically equivalent to that used in 
Ref. In our model the electron-magnon coupling in 
the Hamiltonian Ht (|23|I comes from the exchange inter- 
action of the itinerant electrons mediated by the Coulomb 
potential and is characterized by the matrix element Vq 
Q which depends on the magnon wave-vector q and, 
therefore, on its energy. 



IV. ELASTIC AND INELASTIC 
CONTRIBUTIONS TO THE TUNNELING 
CURRENT 



To evaluate the statistical average of the tunneling cur- 
rent one can use standard first order perturbation theory 
with respect to the transfer Hamiltonian H23|l in the same 
way as in Ref. 0. The voltage drop \e\V between the sys- 
tems is generated by the difference in the Fermi energies 
in the left Ep and the right Ep — \e\V ferromagnets. We 
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will first discuss briefly the elastic contribution Igi : 



lel = 


27Te^VT^ J PMM+P,nrn, P, 
h \ PmM + PAIm, AP, 


(24) 


Pmm 


= PlM {Ef , ^)PrM {Ef , k) , 
k 


(25) 


Pmni 


= 2_^Plm{EF,k.)Prni{EF,iq, 
k 


(26) 


PmM 


= PlmiEp, 'k)prM{EF, k), 
k 


(27) 


PMm 


k 


(28) 



Here pr,iM{EF,ii-) and Pr.hniEp, k) are the majority (M) 
and minority (m) local electron spectral densities related 
to the retarded (R) and advanced (A) electron Green 
functions at the boundaries of the ferromagnets: 

g:^,iEF,k)~g:^,{EF,k) 

Pr,i {Ef , k) = — . (29) 

They can be taken at the Fermi energy for \eV\ ^ Ep- 
As in our case the parallel wave-vector k is conserved 
upon the tunneling (coherent tunneling) , the current l|24|) 



is proportional to the trace of the product of two spectral 
densities. For parallel (P) magnetizations, it is propor- 
tional to pm m + Pmm since the tunneling occurs indepen- 
dently between the majority and minority bands whereas 
in a junction with antiparallel (AP) magnetizations, car- 
riers tunnel from majority to minority bands (and vice 
versa) and hence I^i oc pmM + PMm- The degree of spin- 
polarization V = [I^i — of the elastic current is 
given by 



p _ PMM + Pmm — PmM — PMm ^ ^ ^g^^ 
Pmm + Pmm 



For incoherent tunneling where the parallel momentum 
is not conserved, V would be expressed in terms of 
the local densities of the states J2kPr,iM{EF,^) and 
^i^Pr,imiEF,^) rather than the momentum convolu- 
tions of the spectral densities (|25|l - l|28|) . At the end of 
this paper we give the expressions for the traces H25|l - H28() 
in terms of the band electron parameters of the ferromag- 
nets (see, Eqs. and ^^). 

As to the inelastic current 1^, let us first discuss the 
case of the antiparallel alignment of the magnetic mo- 
ments where one finds 



q ^ ' 

X [{pMM^riuJ, q) -I- pmmfliiuj^q)) {n{e)[l ~ n{e + \e\V - uj)][l + N{uj)] - [1 - n{e)]n{e + \e\V - Oj)N{iu)} + 
+ {pMMni{uj,q) + Prr,mnriuj,q)) {n{e)[l ~ n{e + \e\V + uj)]N [uj) - [1 - n{e)]n{e + \e\V + uj)[l + N (uj)]}] . 



Here we assume that the majority electrons in the left 
and the right systems are spin-up (t) and spin-down (|) 
ones, respectively (Fig. 3), and the magnon spectral den- 
sities are expressed in terms of the advanced and retarded 
magnon Green functions as 



2m 



5{^-u:;'). (32) 



In general, the combinations of the electron n{e) and 
magnon N(uj) occupation numbers in Eq. describe 
four magnon-emission processes and four counter pro- 
cesses involving an absorption of a magnon, which is only 
possible at finite temperatures when N{lo) ^ 0. The 
emission processes are shown schematically in Fig. 3 for 
T = 0. 

The processes in Fig. 3a and 3b correspond to the first 
term in the square brackets in Eq. H31I) which deter- 
mines the current at positive voltages where the Fermi 
energy in the left system is higher than that in the right 
one. The process in Fig. 3a has been already discussed 
in the introduction for the case of half-metallic ferromag- 



nets (Fig. 1). As both initial and final electron states 
belong to the majority bands, the corresponding contri- 
bution to the current is proportional to pmm (I25|) and 
the magnon spectral density at the right side of the junc- 
tion Vlr (|32|l . Note that the similar interaction of the 
minority electrons cannot give rise to the magnon emis- 
sion because it would result in the magnetic moment of 
the right system bigger than that in the ground state. 
However the minority magnon-assisted transport (pro- 
portional to Pmm (EEJ) can be realised in a way shown in 
Fig. 3b: In the course of the exchange scattering, a hot 
minority (spin-down) electron from the left side excites 
a majority (spin-up) electron from the same side above 
the Fermi level of the right system enabling the latter to 
occupy an empty state in the spin-up (minority) conduc- 
tion band in the right ferromagnet. At the same time, a 
spin-wave excitation of the majority (spin-up) Fermi sea 
is created on the left side of the junction, and the corre- 
sponding contribution to the current H31|l is proportional 
to the magnon spectral density at the left boundary Q.i 
The processes in Fig. 3c and 3d are generated in 
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(a) 



(b) 
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^-le+lelV+ODq 
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FIG. 3; Exchange-induced spin- flip tunneling between ferro- 
magnets of arbitrary spin-polarization for antiparallel config- 
uration at zero temperature: (a) and (b) - assisted tunneling 
between the majority and minority bands, respectively, for 
positive bias voltage V > 0, (c) and (d) - same for a nega- 
tively biased junction V < 0. 



27r|5,e|T2 



^ \ Ah 



dedoj X 



(33) 



x{flr{uj, q) -I- Qi{u}, q)) x 

X [pmM {n{e)[l~n{e+\e\V~Lu)][l + N{Lu)]- 
-[1 - n{e)]n{e + \e\V ~ uj)N{uj)\ + 

{n(e)[l - n(e + \e\V + w)]iV(w)- 
-[1 - n(e)]n(e + \e\V + + N{uj)\\\ . 

Unlike the antiparallel alignment (see, Eq. in 
this case the magnon-assisted transport is due to the ex- 
change between a minority electron and a majority one. 
The latter can be either from the opposite or the same 
side of the junction, which explains why the current l|33() 
contains the sum f2r(w, q) -t-r2/(a;, q) of the magnon spec- 
tral densities and is proportional to the trace in the mo- 
mentum space of the product of the minority and the 
majority electron spectral densities H27fl and H28|l . 

In what follows we present the inelastic tunneling spec- 
tra for a simpler case of identical ferromagnets. When 
evaluating the second derivative of the tunneling current, 
one finds that the equilibrium magnon occupation num- 
bers A''(a;) drop out of the expression for (Pl/dV'^. Only 
the non-equilibrium processes of magnon-emission con- 
tribute to the tunneling spectra: 



'2 tAP,P 



dV^ 



T.V^Q 



V2(0) 



de X 



xn'(e) [n'{e + \e\V - Wq) - n'{e + \e\V + lu^)] , (34) 

where the second derivative in the limit of the zero tem- 
perature and voltage determines the amplitude of the 
response: 



2 tAP 



dV^ 



2Tr\S^e^\n 



(35) 



T,V^O 



Ah J 



MM 



d^I^/dV^ 



\T,V^O 



Id^I^P/dV^ 



(36) 



T.V^O 



negatively biased junctions and described by the second 
term in the square brackets in Eq. H31|) . Although they 
look similar to those shown in Fig. 3a and 3b, in general 
there is no symmetry because the magnon spectral densi- 
ties at the left and right boundaries may not be identical: 

For the parallel alignment of the magnetic moments, 
when the majority electrons in both ferromagnets are 
spin-up ones, one can obtain the following expression for 
the inelastic current: 



with fl = A/AttD being the density of states of the sur- 
face magnons. According to Eqs. (|34(l and H36(l . the 
ratio of the amplitudes of the response is controlled by 
the degree of spin-polarization of the elastic current H30|l . 
which can be used for an independent measurement of the 
spin-polarization of the tunneling current. At the same 
time, for identical ferromagnets the shape of the inelas- 
tic tunneling spectrum does not depend on the relative 
alignment of their magnetic moments and the degree of 
spin-polarization (see, e.g. Fig. 2 for half-metallic ferro- 
magnets) . 
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In Fig. 2 both excitation energy \e\V and fc^T are nor- 
malized by the energy of a surface magnon with the wave- 
length equal to the screening radius Dr^^. For T — 
the second derivative is discontinuous at = (curve 
A) recovering the zero-bias anomaly due to the emission 
of surface magnons studied theoretically in Ref. 6. As 
\e\V increases, the wave-length of the excited magnon, 
^ {\eV\/Dy^^ becomes shorter and, when it reaches 
rs, the electron-magnon coupling in Eq. H34() becomes 
strongly energy-dependent leading to a l/|e|V^ decrease 
in the tunneling spectrum. Finite temperatures (curves 
B-E) result in the smearing of the zero-bias anomaly 
due to the counter spin-flip processes (Fig. lb) which 
lead to a finite-rate increase in the response at small 
\e\V and hence to the formation of two antisymmet- 
ric pe aks. In agr eement with the experimental data of 
Refs .13 181121131141 at relatively low temperatures (fc^T < 
Drg^, curves B and C) the peaks are sharp. As the 
temperature increases, they broaden and shift towards 
higher excitation energies (curves D and E). At large \e\V 
all the curves merge showing a temperature-independent 
behaviour, also clearly seen in the experiments^'^. For 
the screening radius rs > a the characteristic magnon 
energy Dr^^ is much smaller than that related to the 
Curie temperature of the ferromagnets ksTc- 

To compare our results with the experimental data 
of Refs. .3.8 12 13-141 we can estimate the voltage cor- 
responding to the peak positions at low temperatures as 
Vp ~ Drg^/|e|. For ~ a ~ lOA and the spin stiffness 
in transition metals D k, 300 — 500 meV x one obtains 
- 3 - 5 mV. In Refs. [H [H and H the peaks were 
observed at 2mV, 12 mV and 17 mV, respectively. At 
the same time, the Curie temperature of the ferromag- 
nets corresponds to the voltage of order of 100 mV. The 
relation (|35|l between the spin-polarization of the current 
and the peak intensities for P and AP configurations is 
also in good agreement with experimental data^^^. 

It remains to calculate the convolutions of the local 
electron spectral densities pmm, Pmm and pMm 
for the isolated right and left systems which are described 
by the boundary conditions H14(l where A means the elec- 
tron penetration length into an infinitely thick barrier 
(a oo) of the finite height U. We will do it for the 
case of identical ferromagnets where it is enough to cal- 
culate the local spectral densities at the boundary of one 
of the electrodes, say, the right one x = a/2. Using the 
expression for the local electron spectral density in terms 
of the advanced and retarded Green functions (|29|l . one 
can write 

PrM,m{EF, k) = J^k^ (g^ (f ) 5{Ef ~ EM,m{k^,k)) = 

^^^EkA9xcf>k^ S{Ef ~ EmM^xM, (37) 

where we introduce the eigenstates (x) in the direc- 
tion perpendicular to the boundary and take into ac- 
count that they must satisfy the boundary condition 
•^fe^ (f ) ~ ^dx4'kj; (§) l(T^ in order to ensure vanish- 
ing of the particle current. For an infinitely high barrier 



(A ^ 0), this boundary condition becomes a "hard wall" 
one and hence 4>k^{x) = (2/ Ly^'^ sink^ix — a/2) with L 
being the length of the system. For a high enough bar- 
rier, we can still use these eigenstates in the second line 
in equation (|37|l since the derivative of sin kx{x — a/2) at 

the boundary is finite. £'Af.m(fcx,k) = ^' ^''■zm'^ ^ ~^ T 
the electron spectrum in the Stoncr model with A mean- 
ing the exchange-induced spin-splitting. Calculating the 
integral in Eq. l|T7|l . one finds 



PrM,m{EF,k) — 



2m{EF ± A/2) 



{kl,,,,-kY/'Q{kM,m-\k\), 

(38) 



where Q{x) is a step- function. The calculation of the con- 
volutions of the local electron spectral densities 
is now straightforward: 



2nm'^A 



Pmm, mm 
PMri 



27:7X1^ A 



-{\kM,mY 



(Afc„ 



(39) 
(40) 



X 



7^/2(7-1-1) /7-I 



arccosh 



7-1-1 
>-l 



7 = kli/kl > 1. 



The same trick with the eigenstates can be used to cal- 
culate the total spin Sz l|18|l of the itinerant electrons 
penetrating into the barrier. It is given at zero tempera- 
ture by 



\^A{k\,i — fc^) 
30(27r)2 ■ 



(41) 



Finally, a strongly nonlinear behaviour of the I — V 
curves at small bias voltages has been recently observed 
in ferromagnet/semiconductor/ferromagnet tunnel struc- 
tures^i. In order to apply our approach to such sys- 
tems one needs to take into account both the energy 
and momentum dependence of the electron penetration 
length into the barrier A since the height of the barrier 
in Ref. |2lj is much smaller than the Fermi energy in the 
ferromagnets unlike the case of oxide junctionsSiiSiiiii. 
We also did not consider interaction processes in which 
a hot electron tunnels elastically through the barrier and 
then emits a magnon in the course of the exchange scat- 
tering in the leadsiP. Since the screening of the Coulomb 
interaction in the leads is much stronger than inside the 
barrier, the discussed mechanism of the nonmonotonic 
behaviour of Cp'I(y)/dV^ may not work in the bulk of 
the ferromagnets. Nevertheless, the proposed theory can 
be extended to account for spin-relaxation in the leads 
through the bias-energy dependence of the local electron 
spectral densities (|29|l . 
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